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We discuss the structure of the nuclear force using a lagrangian derived from hadroniza- 
tion of a chiral quark and diquark model. A generalized trace log formula including meson 
and nucleon fields is expanded to the order in which relevant terms emerge. It is shown that 
the nuclear force is composed of long and medium range parts of chiral meson exchanges 
and short range parts of quark-diquark exchanges. The ranges of the scalar and vector in- 
teractions coincide well with those of sigma (<r) and omega (w) meson exchanges if the size 
of the nucleon core of a quark-diquark bound state is adjusted appropriately. 
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§1. Introduction 



A microscopic description for the meson-nucleon interaction is one of central 
issues in the hadron- nuclear physics. Although lattice calculations seem to reach the 
level of quantitative description of hadron properties, interpretation based on effec- 
tive models is still important in order to understand underlying physics mechanism. 
Symmetry is primarily a powerful tool when one attempts to understand physical 
meaning of some processes. In hadron physics, flavor and chiral symmetries are the 
two important symmetries which one should incorporate. In particular, inclusion 
Qh! of the pion degrees of freedom associated with spontaneous breaking of chiral sym- 

metry is a crucially important ingredient. In addition, for baryon dynamics, their 
internal structure is also important, as various form factors imply. These aspects 
can be intuitively understood by the fact that pion Compton wavelength is larger 
than nucleon size and that nucleon size is larger than nucleon Compton wavelength: 
1/m* » {r 2 N ) 1 ' 2 » l/M N . 

In a recent publication one method was proposed, which can deal with the above 
two aspects of chiral symmetry and internal structure of hadrons using an extended 
model of the Nambu-Jona-Lasinio modeP including not only the quark-antiquark 
but also quark-diquark correlations. 2 ' It was shown that the method works well 
for both mesons and baryons. A microscopic lagrangian was path-integrated to 
generate an effective lagrangian for mesons and baryons with maintaining important 
symmetries such as the gauge and chiral symmetries. Hadron structure was then 
described by their constituents; a quark and an antiquark for mesons and a quark and 
a diquark for baryons. It was also pointed out that the resulting effective lagrangian 
contains various interactions among hadrons, such as meson-meson, meson-baryon 
and baryon-baryon interactions. 

In this paper, we investigate, among many interactions, nucleon-nucleon (NN) 
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interaction*) derived from the previous framework.™ The NN interaction below the 
meson production threshold is phenomenologically well known as derived from the 
phase shift analysis, but its microscopic derivation is still needed. While the long 
range part is well described by meson exchange picture, there are several different 
approaches for the description of the short range part, as including meson exchanges 
and quark exchanges EJSEliElJ3J3 i n our model, all components of the NN force is 
contained in the effective lagrangian which is written in a concise form of trace-log 
type. Then we will see that the NN interaction is composed of the long and medium 
range parts as described by meson exchanges, and short range part of quark-diquark 
exchanges. The latter is then shown to contain various components of the NN 
interaction such as central, spin-orbit and tensor types with non-locality. 

As explained in the next section, there are two types of diquarks; one is scalar, 
isoscalar diquark and the other is the axial-vector, isovector one. In the present 
paper, we perform analytic calculations including the both diquarks and discuss 
general structure of the NN interaction. However, for numerical computation, we 
consider only scalar diquark due to rather complicated structure when the axial- 
vector diquark is included. We then study interaction ranges in more detail than 
interaction strengths, since we expect that the former would be less sensitive to the 
kind of the diquark included. 

We organize this paper as follows. In section 2, we briefly show the derivation of 
the trace-log formula in the path-integral hadronization of the NJL model including 
quark-diquark correlations. In section 3, terms containing the NN interaction are 
investigated in detail, where general structure of the NN amplitude is presented. We 
show a sample numerical calculation only for the case containing a scalar diquark. 
The present study of the NN interaction is not quantitatively complete but will be 
useful in showing some important aspects of the nuclear force, in particular that the 
range of the short range interaction is related to the intrinsic size of the nucleon. 
The final section is devoted to summary. 

§2. Effective lagrangian for mesons and nucleons 

We briefly review the method to derive an effective lagrangian for mesons and 
nucleons from a quark and diquark model of chiral symmetry following the previous 
work by Abu-Raddad et alP^ Let us start from the NJL Lagrangian 

G 

£njl = q(ifi ~ m )q + — [{qqf + (qi^rq) 2 ] ■ (2-1) 

Here q is the current quark field, f the isospin (flavor) Pauli matrices, G a dimensional 
coupling constant, and uiq the current quark mass. In this paper we set mo = for 
simplicity. As usual, the NJL lagrangian is bosonized by introducing collective meson 
fields as auxiliary fields in the path-integral method!*^'^'^ At an intermediate step, 



*' In what follows we use "nucleon" rather than "baryon", as we consider flavor SU(2) for the 
nucleon sector. 
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we find the following lagrangian: 



1 



C' gan = q (ip - (a + i 75 r • if)) q - ^ (a 2 + if 2 ) . (2-2) 

Here a and 7r are scalar-isoscalar sigma and pseudoscalar-isovector pion fields as 
generated from a ~ qq and n ~ iqr^q, respectively. For our purpose, it is convenient 
to work in the non-linear basis rather than linear oneP^J 1 ^ It is achieved by the 
chiral rotation from the current (q) to constituent (x) quark fields: 

/ . ■ -> ->\ 1/2 

X = £5<?, &=( j J , (2-3) 

where f 2 = a 2 + tt 2 . Thus we find 



where 



r' vrrr =xW-f-t-t)x--^f 2 , (2-4) 



p 2 

are the vector and axial-vector currents written in terms of the chiral field 

&=(^P)" 2 (2-6) 

The lagrangian (|2-4|) describes not only the free kinetic term of the quark, but also 
quark-meson interactions such as Yukawa type, Weinberg- Tomozawa type and etc. 

In a model we consider here, we introduce diquarks and their interaction terms 
with quarks. We assume local interactions between quark-diquark pairs to generate 
the nucleon field. As inspired from the method to construct local nucleon fields, it is 
sufficient to consider two diquarks; one is scalar, isoscalar one, D, and the other is 
an axial-vector, isovector one, D^!^ Hence, our microscopic lagrangian for quarks, 
diquarks and mesons is given by 

C = m -f-t-t)x- 7^f + D\d 2 + M 2 )D + " [{d 2 + M\) 9ilv - d^d v ] & 
+ G ( sin 9 xi^l 5 t " 

fi + cosdxD ] ^ (sm9 D u -t^1 5 X + cos6 D X ) • (2-7) 



In the last term G is another coupling constant for the quark-diquark interaction 
and the angle 9 controls the mixing ratio of the scalar and axial- vector diquarks in 
the nucleon wave function. 

One would have constructed a nucleon field in the linear basis. However, this 
makes chiral transformation properties complicated. In fact, if we replace x by 
q, then more than two terms are necessary to maintain chiral symmetry for the 
interaction term in (|2-7|) . In contrast, in the non- linear basis, the transformation 
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property becomes simple; both quark and nucleoli fields are isospin multiplets which 
are subject to the non- linear transformation 

X{x) -> h(x) X (x) , B(x) -» h(x)B(x) . (2-8) 

Here h is a non-linear function of the chiral transformations and the chiral field at a 
point x.OJEI 

Now the hadronization procedure, the elimination of the quark and diquark fields 
in (|2-7|) . is straightforward. The final result is written in a compact form aa^ 

^ e ff = -7^/ 2 " itrln(^-/-^-^) - IBB + i tr ln(l - □) . (2-9) 

Here trace is over space-time, color, flavor, and Lorentz indices, and the operator □ 
is defined by 

5), (2-10) 

where 

A 1 *' uj = sm 2 9 B 7p7 5 r k A pk '^ S r j 7 V B , (2- 11a) 

5 = cos 2 6> B AS B , (2-llb) 

= S m9cos9B A S T j 7 V 5 , (2-llc) 

(^j,)^ = sin^ cos S i^'^ 7p7 5 r fc SB. (2- lid) 

Though the effective meson-nucleon lagrangian (|2-9|) looks simple, it contains 
many important physics ingredients when the trace-log terms are expanded: 

• It generates a meson lagrangian in a chirally symmetric manner. Up to the 
fourth power in the meson fields, it produces precisely the lagrangian of the 
linear sigma model with realization of spontaneous breaking of chiral symmetry. 

• From the second trace-log term, a nucleon effective lagrangian is derived. In 
the previous paper, kinetic term of the nucleon was investigated and nucleon 
mass was computed at the one-loop levelP^ 

• In the nucleon effective lagrangian, meson-nucleon couplings appear through the 
diagrams as shown in Fig. ^ Their strengths and form factors are computed by 
the underlying quark-diquark dynamics. Using these vertices, meson-exchange 
interactions are constructed. 

• There are diagrams which contain many nucleon fields. For instance, NN in- 
teractions are expressed by one-loop diagram as shown in Fig. |21 This term 
describes a short range part of the NN interaction. In this paper, we focus our 
attention mostly on the NN interaction derived from the one-loop diagrams. 

§3. The structure of the NN interaction 



As we have explained, the NN interaction in the quark-diquark model consists of 
two parts: one is the meson exchange terms for long and medium range interactions 
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Fig. 1. Meson- nucleon Yukawa vertex. 




P'l 




Fig. 2. A loop diagram for the NN interaction (left) and the equivalent quark exchange 



and the other is the one described by the quark-diquark one-loop diagram shown 
in Fig. [21 for the short range interactions. The latter diagram is equivalent, when 
the flow of nucleon lines are arranged appropriately, to quark or diquark exchange 
diagrams as shown in Fig. |3 In this quark exchange process, the nucleon goes 
into the quark-diquark pair through the interaction (blobs in the diagram). In this 
sense, this process differs from the norm kernel of quark exchange in a quark cluster 



3.1. One pion exchange interaction 

Let us first examine the long rage part which is dominated by the one pion 
exchange. In the present model setup, we start from the Nambu-Jona-Lasinio model 
which includes the pion and sigma mesons. In the nuclear force, the scalar-isoscalar 
interaction also emerges from quark exchanges at the same energy scale as the sigma 
meson exchanges. This is a general remark when we consider the medium and short 
range parts of the nuclear force; there are two contributions from boson exchanges 
and quark exchanges, a natural consequence when nucleon consists of quark and 
diquark core surrounded by meson clouds. 

Now the pion exchange process in the present model is supplemented by the 
form factor due to the quark-diquark structure of the nucleon as shown in Fig. ^ In 
the momentum space, the one pion exchange interaction is then written as, in the 
non-relativistic approximation, 



diagram. 



mode lP,EJ 



VoPEp(^) 





2M N 



q 2 + m 2 . 



2M N 



where g-wNNiq 2 ) is the irNN form factor, Mn the mass of the nucleon and are 
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the spin matrices for the nucleoli 1 and 2. Due to chiral symmetry, it is related to 
the axial form factor through the Goldberger-Treiman relation^ 

M N 

9-kNN = —j—9A ■ (3-2) 

Jtt 

In the present paper, we include only the scalar diquark due to rather com- 
plicated structure of the quark-diquark loop integral when the axial-vector diquark 
is included. A full calculation with both scalar and axial-vector diquarks included 
requires careful treatment of the regularization procedure. We will postpone this 
complete study in future works. It is, however, not very difficult to estimate the ef- 
fect of the axial-vector diquark for one-body matrix elements; it mostly contributes 
to the absolute values of spin-isospin (a ■ r) matrix elements, which are axial- vector 
coupling constant gA (and hence the irNN coupling constant g w NN) arid isovec- 
tor magnetic moments. For quantities which reflect the size of the quark-diquark 
bound state such as q 2 dependence of form factors are not very sensitive to the kind 
of diquarks, since they are dominated by the extension of the lighter particle, the 
quark. 

To start with, numerical calculations are performed using the parameters as 
shown in Table [IJ These were used in the previous work for the study of a single 
nucleon. 2 ' The divergent integrals are regularized by the Pauli-Villars method with 
the cut-off mass .Apy- Then we perform calculations by varying the coupling 
constant for the quark-diquark interaction G. The strength of G controls the binding 
or size of the nucleon as shown in Table 1. Here the size is defined as a sum of the 
quark and diquark distributions: 

(r 2 ) = (r 2 g ) + (r 2 D ). (3-3) 

Note that this is different from the charge distribution where each term is weighted 
by the corresponding charge. We will see how physical quantities such as form factors 
and short range interactions are related to the extension of the quark-diquark wave 
functions. 



Table I. Model parameters taken from Ref. 



Mn Vtlq 


M D 


Apv 


G 


0.94 GeV 0.39 GeV 


0.60 GeV 


0.63 GeV 


271.0 GeV" 1 



In Fig. |3J we plot the ttNN form factor as a function of the square momentum 
transfer q 2 for different binding energies. We did not attempt to reproduce the 
experimental value of the nucleon mass, since our interest here is how the form factor 
is related to the extension of the nucleon wave function due to the quark-diquark 
structure. The result is the same as the previous work when the parameters shown 
in Tableware used. As G is increased such that the size of the nucleon is reduced, 
the g-dependence of the form factor becomes weak. We have extracted the cut-off 
parameter of the form factor when fitted by a monopole function, A 2 /(A 2 + q 2 ), 
and the range of the axial form factor, (r 2 ) 1 ^ 2 . The result is shown in the second 
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Table II. Basic nucleon properties when G is varied. Other model parameters m q , Md and A are 
fixed. 



G(GeV _1 ) 


B.E.(MeV) 


Miv(MeV) 


(rT^(fm) 


156.4 


10 


980 


1.64 


271.0 


50 


940 


0.89 


445.9 


140 


850 


0.60 


529.4 


190 


800 


0.54 


609.2 


240 


750 


0.51 


687.4 


290 


700 


0.48 



and third column of Table II I II When the previous parameter set is used such that 
the nucleon (electromagnetic) size about 0.8 fm is reproduced by the quark-diquark 
wave function^ the form factor becomes rather soft, A ~ 500 MeV. However, if 
we consider the fact that meson clouds should also contributes to the nucleon size 
and the core part should be smaller than the observed size, then the quark-diquark 
extension could be smaller. If we take, for instance, the mean square value of the 
quark-diquark extension about (0.6 fm) 2 , then the cut-off size of about 800 GeV 
is obtained. This value is smaller as compared with those often used in the NN 
forceP^J However, the cut-off parameters in the NN force depends on the modeling 
of the meson exchanges. In fact, the above value is consistent with that extracted 
from the electroproduction of the pionP^ 



Table III. The cut-off parameter and the corresponding range of the uNN form factor. 



Mjv(MeV) 


A(GeY) 


(^)^(fm- 1 ) 


980 


0.29 


1.67 


940 


0.54 


0.90 


850 


0.78 


0.62 


800 


0.85 


0.57 


750 


0.90 


0.53 


700 


0.93 


0.52 



3.2. Short range interaction 

Let us turn to the short range interaction described by the quark (or diquark) 
exchanges, as shown in Fig. [3 Using the interaction vertices given in Eq. (|2-10|) . it 
is straightforward to compute the amplitude for the quark-diquark loop: 

M NN = -iN c Z 2 

x [_#k B(t + m q )BB(t 2 -j>' 2 + t + m q )B 

J (2vr) 4 [{ Pl - k) 2 - M 2 ] [k 2 - m 2 } [(p' 2 - k) 2 - M 2 } [{p 2 - p' 2 + k) 2 - m 2 } ' 1 ' 

where N c is the number of colors and Z is the wave function normalization of the 
nucleon fieldP' In this paper, we evaluate the integral in the center-of-mass system 
for the elastic scattering of the two nucleon: 

Pi = (Ep, p), P2 = {E$, -p) 

Pi = (■%'> P'), P2 = i E p', -p') , \p\ = \p'\ ■ 
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0.0 1 1 1 

0.0 0.5 1.0 

q 2 [GeV 2 ] 

Fig. 3. Pion-nucleon form factor g^NN- Three lines are for M/v=0.98 (loosely bound), 0.85, 
0.70 (tightly bound) GeV from bottom to top. 

In principle, from this amplitude we should calculate observables such as phaseshifts, 
cross sections for physical two-nucleon channels by properly taking into account 
anti-symmetrization. Since in the present calculation we include only the scalar 
diquark, we do not expect that the comparison of such quantities with data will 
make significant sense. Rather, in the following we study some basic properties of 
the amplitude itself, mostly the interaction ranges extracted from Eq. (j3-4j) . 
To do so, we write the amplitude (|3-4|) in a schematic manner as 

M NN = F s (P,q)(BB) 2 + F v (P,q)(B lfl B) 2 + -- - , (3-5) 

where 

P = p' + p q=p'—p (3 - 6) 

In Eq. (j3-5[l . we identify Fs and Fy with the scalar and vector interactions, respec- 
tively. Dots of Eq. (|3-5|) then contain terms involving external momenta pi such as 
Br(pi)B, where r(pi) is a 4 x 4 matrix involving p, L . The latter are usually not con- 
sidered in the local potential but here they exist due to the non-local nature of the 
amplitude (|3-4[1 . Furthermore, we note that the bilinear forms of the nucleon fields 
in Eq. (|3-5|) are taken by the pair of B fields in Fig. 21 carrying the momentum p\ and 
p' 2 , and the other carrying the momentum p 2 and p^. This is a consequence of the 
quark-diquark loop integral. Before going to details, it is pointed out that the scalar 
interaction is attractive, while the vector interaction is repulsive. This is a general 
feature of the interaction between two nucleons (fermions) and is independent of the 
details of a model. 

In general the amplitude is highly non-local as the symmetric loop diagram 
implies. The loop diagram may be interpreted as either a quark exchange or diquark 
exchange diagram. It is then natural that the heavier diquark exchange is of shorter 
range nature. By identifying the momentum transfer exchanged by a quark pair as 
q and that of the diquark pair by P, it is expected that P dependence is milder than 



Structure of the nuclear force 



9 



q dependence. We show this explicitly in Fig.0J Therefore, we expand the functions 
Fg v m powers of P and evaluate their coefficients as functions of q: 



F i (P,q)=F i (P,q)\ 



+ 



(3-7) 



where i stands for S or V. The resulting q dependent functions might be interpreted 
as a Fourier transform of a local potential as a function of the relative coordinate 
f = x\ — %2- Then the general structure of the NN potential can be studied by 
performing non-relativistic reduction of the amplitude, Eqs. (|3-5|) or (|3-7|) , It contains 
central, spin-orbit, tensor components as well as non-local terms as proportional to 
P. 



0.6 



0.0 





F s (P 2 ,0) 




s F s (0,0) 


F s (0,q 2 ) 




' f s (0,0) 





0.0 



0.5 

q or P 2 [GeV 2 



L.O 



0.8 



0.4 





F v (P 2 ,0) 
F v (0,0) 


F v (0,q 2 ) 




Fy (.0,0) 



0.5 

q 2 or P 2 [GeV 2 ] 



Fig. 4. Normalized scalar and vector form factors as function of q 2 = \q 2 \ with fixed P = 
(solid line) and those of P 2 = \P 2 \ with fixed q = (dashed line). Calculations are 
performed using the parameter set of Table 1. 



Now we discuss the scalar- isoscalar functions F$ and Fy in Eq. (|3-5j) . We can not 
write the resulting q dependence in a closed form, but have to resort to numerical 
computations. We have performed the Feynman integral numerically for several 
different G parameters which yields different binding energies, or the size of the 
quark-diquark distribution. Although the one-loop integral of Fig. [2 converges for a 
scalar diquark, we keep the counter term of the Pauli-Villars regularization. Since the 
present theory is a cut-off theory with the relatively small cut-off mass Apy = 0.63 
GeV, the counter term plays a significant role in order to produce numbers. 

First we see the strength of the interaction, by extracting the Yukawa coupling 
constant of the vector and scalar type from the amplitude. The results are shown in 
Fig. [3] as functions of the size of the nucleon (r 2 ). These numbers may be compared 
with the empirical values gs ~ 10 and gy ~ The present results are strongly 

dependent on (r 2 ). When only scalar diquark is included the scalar interaction 
becomes much stronger that the strength of the vector interaction. Phenomenologi- 
cally, the vector (omega meson) coupling is slightly stronger than the scalar (sigma 
meson) coupling. We expect that the present result will be changed significantly 
when axial-vector diquark is included. 

Let us turn to the discussion of the interaction ranges which is expected to 
be more reliably studied than the strength even when only the scalar diquark is 
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0.4 0.6 0.8 1.0 1.2 1.4 
<i- 2 >" 2 [fm] 



Fig. 5. Scalar and vector coupling constants as functions of the nucleon size (r 2 ) 1 / 2 . 



included. In Fig. H3 we have shown the (/-dependence of the zeroth order coefficients 
of (|3-7j) for three different sizes of the nucleon. It is obvious from Fig. H3 as the size 
of the nucleon becomes smaller the interaction range becomes shorter range. More 
quantitatively, we define the interaction range Ri by 



-6 



1 dFi 



Fi(q 2 )dq 2 



(3-8) 



which is related to the mass parameter of the interaction range by rrii = \/6/i2j. It 
is interesting to see that the mass (and hence range) parameters of the interaction 
are approximately proportional to the size of the nucleon (r 2 ) 1 / 2 as shown in Fig. [71 





0.5 
q~ [GeV 2 



1.0 



0.0 



0.5 

q [GeV 2 ] 



Fig. 6. Normalized scalar and vector form factors Fs,v(q 2 , P 2 = 0)/Fs,v{q 2 = 0, P 2 = 0). 
Each lines are for 1/^=0.98, 0.85, 0.70 GeV from bottom to top. 



When the previous parameter set is used, the interaction mass parameter is 
about 550 MeV and 600 MeV for vector and scalar interactions, respectively. How- 
ever, if we use a parameter set that produces the nucleon size about 0.6 fm, then the 
two masses are ms ~ 600 MeV and my ~ 800 MeV. These numbers are remarkably 
close to the mass of the sigma and omega mesons. The nucleon size of 0.6 fm should 
be reasonable as the size of the nucleon core. Meson clouds around it contributes to 
a substantial amount to the size of the nucleonP^ 

In the present analysis, apart from the absolute value of the interaction strengths, 
the interaction ranges for the scalar and vector interactions have been produced 
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3.0 
2.5 
2.0 
1.5 
1,0 
0.5 
0.0 



1.0 

[fm] 




Fig. 7. The interaction range Ri (left panel) and the corresponding mass parameter m s 
(right panel) as functions of the nucleon size (distribution of a quark and a diquark.) 



appropriately by the quark exchange picture. One question concerns the small (but 
non-negligible) difference between the ranges in the scalar and vector channels, which 
is consistent with the empirical fact. It is crudely understood by the dimensionality 
of the loop integral. As seen from Eq. (J3-4j) the integrand for the vector interaction 
has higher power in momentum variable than the scalar one. Because of this, the 
vector part reflects shorter-distant dynamics and produces shorter interaction range. 

§4. Summary 

In this paper we have studied the nucleon-nucleon (NN) interaction using a mi- 
croscopic theory of quarks and diquarks. Nucleons are described as quark-diquark 
bound states, which was shown to explain reasonably well their static properties.^ 
The quark and diquark degrees of freedom were integrated out in the path-integral 
method and an effective lagrangian was derived for mesons and baryons. The result- 
ing trace-log formula contains various meson and nucleon interaction terms, including 
the NN interaction in the short range region expressed by a Feynman integral of a 
quark-diquark loop. Hence the nucleon-nucleon interaction is naturally expressed as 
meson exchanges at long ranges, while quark and diquark exchanges at short ranges. 
The NN interaction contains the scalar and vector type terms at leading orders in 
powers of momentum of the Feynman integral. It turns out that the scalar term is 
attractive, while the vector term repulsive. 

In the present paper our numerical calculations contained only scalar diquarks 
as a first step toward full calculations. Therefore, the magnitude of the interaction 
strengths were not studied; they must be affected by the inclusion of the axial- vector 
diquark. On the other hand the interaction range was better studied in the present 
framework, reflecting the size of the nucleon. Consequently, the range of the scalar 
and vector interactions came out to be quite reasonable; in the mass parameter, 
the scalar mass was about 600 MeV and the vector mass about 800 MeV, when the 
nucleon size was set at around 0.6 fm. 

The present result encourages us to study further baryon properties by extending 
the model by including the axial-vector diquark. We have already started the study 
of such direction. When the axial vector is included, due to the massive vector nature 
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of the propagator the loop integral diverges more strongly than the case of the scalar 
diquark. This complicates numerical study more than the present case. We hope to 
make progress along this line in the future. 
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